Fermions in gravity and the skyrmion backgrounds in six dimensional brane-worlds 
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We construct brane solutions in six dimensional Einstein-Skyrme systems. A class of baby 
skyrmion solutions realize warped compactification of the extra dimensions and gravity localiza- 
tion on the brane for negative bulk cosmological constant. Coupling of the fermions with the brane 
skyrmions successfully lead to the brane localized fermions. The standard representation of the 
gamma matrices is used to obtain massive localized modes as well as the massless one. Nonlin- 
ear nature of the skyrmions brings richer information for the fermions level structure. In terms of 
the level crossing picture, emergence of the massive localized modes as well as the zero mode are 
observed. 
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I. INTRODUCTION 

Theories with extradimensions have been expected to 
solve the hierarchy problem and cosmological constant 
problem. Experimentally unobserved extradimensions 
indicate that the standard model particles and forces are 
confined to a 3-brane [ll, 0, 0] ■ Intensive study has been 
performed for the RS brane model in five space-time di- 
mensions 2, 3] . In this framework, the exponential warp 
factor in the metric can generate a large hierarchy of 
scales. This model, however, requires unstable negative 
tension branes and the fine-tuning between brane ten- 
sions and bulk cosmological constant. 

There is hope that higher dimensional brane models 
than five could evade those problems appeared in five di- 
mensions. In fact brane theories in six dimensions show 
a very distinct feature towards the fine-tuning and nega- 
tive tension brane problems. In Refs. [J, [^], it was shown 
that the brane tension merely produces deficit angles in 
the bulk and hence it can take an arbitrary value without 
affecting the brane geometry. The model is based on the 
spontaneous compactification by the bulk magnetic flux. 
If the compactification manifold is a sphere, two branes 
have to be introduced with equal tensions. If it is a disk, 
no second 3-brane is needed. But still the fine-tuning be- 
tween magnetic flux and the bulk cosmological constant 
can not be avoided although non-static solutions could 
be free of any fine-tuning i&\ . 

Alternatively to the flux compactification in 6 dim., 
the nonlinear sigma model has been used for compactifi- 
cations of the extra space dimensions [3, M, S [l3l ■ ^^ in 
the flux compactification, no second 3-brane is needed if 
the parameters in the sigma model and bulk space-time 
are tuned. 

Warped compactifications are also possible in six 
space-time dimensions in the model of topological ob- 
jects such as defects and solitons. In this context 
strings [iil, Il2j [S [3 [3 were investigated, showing that 
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they can realize localization of gravity. Interestingly, if 
the brane is modeled in such a field theory language, the 
fine-tuning between bulk and brane parameters required 
in the case of delta-like branes turns to a tuning of the 
model parameters [16| . 

The Skyrme model is known to possess soliton solu- 
tions called baby skyrmions in two dimensional space [13, 
flSll . In this paper we therefore consider the warped com- 
pactification of the two dimensional extra space by the 
baby skyrmions. We find that in the 6 dim. Einstein- 
Skyrme systems, static solutions which realize warped 
compactification exist for negative bulk cosmological con- 
stant. Since the solution is regular except at the coni- 
cal singularity, it has only single 3-brane. Thus no fine- 
tuning between brane tensions is required. The Skyrme 
model possess a rich class of stable multi soliton solutions. 
We find various brane solutions by such multi-solitons. 

It should be noted that general considerations in the 
6 dim. brane model with bulk scalar fields suggest that 
the mechanism of regular warped compactification with 
single positive tension brane is not possible [l^. How- 
ever, the model under consideration is restricted to the 
bulk scalar field depending only on the radial coordinate 
in the extra space. The scalar field in the Skyrme model 
depend not only the radial coordinate but also the angu- 
lar coordinate to exhibit nontrivial topological structure, 
which makes possible to realize regular warped compact- 
ification. 

Study of localization of fermions and gauge fields on 
topological defects have been extensively studied with 
co-dimension one [13, [SI, [H, [H and two [13, [H, [H, 
Wn . [28l [2g |. In many years before, particle locahzation 
on a domain wall in higher dimensional space time was 
discussed [sol |31|. The authors suggested the possibility 
of localized massless fermions on the 1 dim. kink back- 
ground in 4-t-l space-time with Yukawa- type coupling 
manner. Localization of chiral fermions on RS scenario is 
in Rcf. 20] . Analysis for the massive fermionic modes was 
done by Ringeval et, at in Ref. 21] . For co-dimension two, 
the localization on higher dimensional generalizations of 
the RS model was studied by the coupling of real scalar 
fields [2^. Many studies have been followed and most of 



them are based on the Abehan Higgs model coupled with 
the chiral fermions. 

Problem of fermion mass hierarchy was discussed in 
Ref . [25l. |26| . |33 . 33 1 within different mechanisms. Es- 
pecially, in Ref. [2^ the hierarchy between the fermionic 
generations are explained in terms of multi- winding num- 
ber solutions of the complex scalar (Higgs) fields. They 
observed three chiral fermionic zero modes on a topolog- 
ical defect with winding number three and finite masses 
appear the mixing of those zero modes. Although any 
brane localization mechanism is absent in their discus- 
sion, the idea is promising. 

In this article, we employ somewhat different set up: 
we consider the localization of the fermions on the 
skyrmion branes and, to treat the massive fermionic 
modes directly, we use the standard representation of 
the higher dimensional gamma matrices instead of the 
chiral one. The fermion localization and the existence 
of the zero modes arc confirmed through the analysis of 
spectral flow of the one particle state. 

This paper is organized as follows. In the next section 
we describe the Einstein-Skyrme system in six dimen- 
sions and derives the coupled equations for the Skyrme 
and gravitational fields. We derive a class of multi- 
winding number solutions. We will show some typical nu- 
merical solutions. Formulation of the fermions in higher 
dimensional curved space time is discussed in Sec. III. 
Coupling of the fermions and the skyrmions is introduced 
in this section. Conclusion and discussion are given in 
Sec.IV. 
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where M,N run over 0, • • • ,5 and n = (0,0,1). 4> = 
{(j)^, (fP' ,4''^) denotes a triplet of scalar real fields with the 
constraint (j) ■ (f> = 1. The F,e,^ are the Skyrme model 
parameters wih the dimension of M^, M~^, M^ respec- 
tively. The first term in ^ is familiar from cr— model. 
The second term is the analogue of the Skyrme fourth 
order term in the usual Skyrme model which works as 
a stabilizer for obtaining the soliton solution. The last 
term is referred to as a potential term which guarantee 
the stability of a baby-skyrmion. 

The solutions of the model would be characterized by 
following topological charge in curved space-time 
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where V^ means the space-time covariant derivative. Let 
us assume that the matter Skyrme fields depend only on 
the extra coordinates and impose the hedgehog ansatz 



(f) = {sin f (r) cos nd , smf{r)smn9, cos/(r)). 



(6) 



The function /(r) which is often called profile function, 
has following boundary condition 



/(O)--(to-I)^, lim /(r) 
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II. CONSTRUCTION OF THE 
BABY-SKYRMION BRANES 

A. Model 

We consider the model of the 6 dim. Einstcin-Skyrme 
system with a bulk cosmological constant coupled to 
fermions. The action comprises 
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Here S'gravity is the six dimensional Einstein-Hilbert ac- 
tion 
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In the parameter k^ — 1/M|, Mq is the six-dimensional 
Planck mass, denoted the fundamental gravity scale, and 
Ab is the bulk cosmological constant. 

For S'brano wc usc the action of baby-Skyrme model [17|, 
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where (tti, n) is arbitrary integer. This ansatz ensures the 
topological charge 



g = n(l-(-l)™)/2. 



(8) 



We consider the maximally symmetric metric with van- 
ishing 4D cosmological constant, 



ds^ = B^{r)r]^^dx''dx'' + dr^ + C^{r)d0^ 



(9) 



where 77^^ is the Minkowski metric with the signature 
(— , -|-, +, +) in our convention and < r < 00 and < 
6 < 27r. This ansatz has been proved to realize warped 
compactification of the extra dimension in models where 
branes are represented by global defects [3J]. 

'S'fcrmion is the actiou of a fermions coupled with the 
skyrmions and the warp factors; that would be described 
in Sccimi 

The general forms of the coupled system of Einstein 
equations and the equation of motion of the Skyrme 
model are 
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where the stress-energy tensor Tmn is given by 
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Inserting Eq. © into Eq. (JH) , one obtains the Lagrangian 
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where we have introduced the dimensionless quantities 
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and 



u = 1 + 



9 • 2 /• 

n sm J 
C2 



(15) 



The prime denotes derivative with respect to the radial 
component y of the two extra space. The Skyrme field 
equation is thus 
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where 
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Within this ansatz, the components of the stress-energy 
tensor (fT2l) becomes 
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The Einstein equations with bulk cosmological constant 
are written down in the following form 

3b' + 6P + 3bc + c' + 0^ = -a{Ab + To{y)) (22) 
6b^ + 4bc^~a{Ab + Tr{y)) (23) 

4b' + 10b'' = -a{Ab + Te{y)) (24) 



where a = n'^F'^ is a dimensionless coupling constant 
and A;, — K^/e^F'^ is a dimensionless bulk cosmological 
constant. Also, we introduce b = B' /B, c = C /C for 
convenience. 



B. Boundary conditions 

At infinity, all components of the energy-momentum 
tensor vanishes and the Einstein equations ([22l) - (j24p are 
then reduced to 

36' + 66^ + 36 c + c' + c^ = -aA;, (25) 

Qb' + Abc^-aAb (26) 

46' + 106^ = -akb . (27) 

The general solution has been obtained in Ref. [ij, [l^ 
which is given by 



b = p ? p — , c = — 

Ae^py + e-^Py 26 2 

where A is an arbitrary constant and 
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Since we are interested in regular solutions with warped 
compactification of the extra-space, the functions B and 
C must converge at infinity. This is achieved only when 
A(, < and ^ = with the solution 
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where Ci and 62 are arbitrary constants. Then, the 
asymptotic form of the metric which realizes warped com- 
pactification is given by 
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The four-dimensional reduced Planck mass Mpi is de- 
rived by the coefficient of the four-dimensional Ricci 
scalar, which can be calculated inserting the metric Q 
into the action 0, 

^ /d^xx/^^i^W = ^ f d^xV^B-'ir)R^^^ 

^ Me. f d^xy/^^R^^^ f drd9B'{r)C{r) 

^^^^«' f drB\r)C{r) f d^x^/^^R^^^ 



where the superscript (4) represents a tensor defined on 
the four-dimensional submanifold. Thus, we find the re- 
lation between Mpi and Afg as 



M^i = 2nMt 



drB'{r)C{r). 
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FIG. 1: Typical results of the profile functions / FIG. 2: Typical results of the stress-energy tensors ro, rr, re 
(straight line), the warp metrices B,C (dashed,dotted (dashed, dotted, and dot-dashed line, respectively) and the 
line, respectively), as a function of y. topological charge density q (straight line). 



The requirement of gravity localization is equivalent to 
the finiteness of the four-dinicnsional Planck mass. For 
the asymptotic solution (pTj) . the localization is attained. 
Let us consider the asymptotic solutions for skyrmions. 
we can write 



f{y)^f + 5f{y). 



(33) 



where for y ^ 1, / ~ 0. The linearized field equations 
are given by 



5f"~^p5f'~ii5f^0. 



(34) 



Assuming that / falls off exponentially, one obtains for 
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where /c is an arbitrary constant. 

Following regularity of the geometry at the center of 
the defect are imposed 
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B'{0) = 0, C(0) = 0, C"(0) == 1 



(36) 



and we can arbitarily fix B{0) = 1. Boundary conditions 
for the warp factors and the profile function at the origin 
are determined by expanding them around the origin. 
For the different topological sectors, the first few terms 
are schematically written down as 

/(y) - -(m - l)7r + /("Ho)y" + 0(y"+i) (37) 
biy)^By + Oiy') (38) 

(7(y)=y + Cy3 + 0(y5) (39) 



Thus one finds that the only /'(O) or /"(O) is the free 
parameter vicinity of the orgin. 



Consider linear combinations of Eas. p2)l - ((24|) . we ob- 
tain 
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Integrating Eqs. (|44| . (|45| from zero to Xc, we get 
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B^Cdx - -{Afio + Air - A^e) . (47) 



(|15)l is the six-dimensional analogue of the relation deter- 
mining the Tolman mass whereas Ea. (|47p is the general- 
ization of the relation giving the angular deficit. Com- 
bining these the following relations are obtained in the 
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These conditions are used for checking the numerical ac- 
curacy of our calculations. 



C. Numerical solutions 

The equations (ITn]),(I221)-(I21]) should be solved numer- 
ically since they are highly nonlinear. The simple tech- 
nique to solve the Einstein-Skyrme equations are the 
shooting method combined with the 4th order Runge- 
Kutta forward integration |35j . However, a unique set 
of boundary conditions at a; = produces 2 distinct so- 
lutions, one of which grows exponentially and the other 
decays exponentially as a; — > cx). This causes instability of 
solutions when the forward integration is performed. In- 
stead, we use a backward integration following Refs. [IJI 
where the 6 dim. vortex-like regular brane solutions were 
constructed. The backward integration method requires 
a set of boundary conditions at infinity. We, however, 
truncate and take the distance x^ax far enough from 
the origin so that the Skyrme profile would fall off be- 
fore it reaches Xmax- The set of boundary conditions at 
X — Xmax produces a unique solution which satisfies the 
boundary conditions at x = and hence it is numeri- 
cally stable. We present our typical numerical results in 
Figs.IHa. 

III. FERMIONS 

A. Basic formalism 

The action of the fermions coupled with the Skyrme 
field in a Yukawa coupling manner can be written as 



*^fermion — / ^ •^*-'ferinion 



(50) 



with 



L 



fcnnion 



= y^ "^{iT^Da - Mt-(P)^ 



(51) 



The six-dimensional gamma matrices F are defined with 
the help the vielbein ef and those of the flat-space 7°, i.e., 
r^ := e^7°. The covariant derivative is defined as 
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S/ Ae% are the spin connection with 
i [7a : 7s] • "^^^ simbol d implies 



generators ct-i^ 

that ^pl)(t) = il^dcj) - (ai/')0. Here A,B = 0, • • • , 5 
are the six-dimensional space-time index and a,b = 
0, • • • ,5 corresponds to the flat tangent six-dimensional 
Minkowski space. The vielbein is defined through qab — 

in Refs.[26|,[27| is simply defined as 



e^eas = ^JsteleR. The definition which was introduced 
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They are the definitions which produce the gamma ma- 
trices parallel to the polar coordinate {r,0). Thus, spe- 
cial care is needed to explore the conserved quantities 
like angular momenta from the corresponding hamilto- 
nian. Non-vanishing components of the corresponding 
spin connection are found to be 

c.f = <5^i3'(r), 4^^-C'{r) (54) 

where ' = ^. The Dirac equation takes the form 
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The corresponding hamiltonian reads 



H = -^^S\^r -f ^ + g) - z7S^^a« + 7°Mf • $. 

(56) 

The equation (|55p can be solved numerically in terms 
of, for example, the shooting method. At present study, 
however, we will treat the problem by somewhat different 
way. According to the Rayleigh-Ritz variational method 
|3a | , the upper bound of the spectrum can be obtained 
from the secular equation; 
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FIG. 3: Fermion number density for the background skyrmion 
with {m,n) — (1,1). The results for the coupling constant 
m = (decoupled) and m — 0.86 are plotted. If fermions 
couple to the skyrmions, only the state Ks = is localized on 
the brane core. The all other states are not observed because 
they are strongly delocalized. 




FIG. 4: Fermion number density for the background skyrmion 
with (m,n) = (1, 1) for the varying the coupling constants. 




where {iy9i} [i = 1, • • • , N) is some complete set of the 
plane- wave spinor basis. For N —i- oo, the spectrum e 
becomes exact. Eg. ([57]) can be solved numerically. For 
simplicity, we are to construct plane- wave basis in the flat 
space-time, i.e.,B ~ 1,B' = 0,C = r, as {(fii}- However, 
no corresponding flat Hamiltonian subject to Eq. (j56p ex- 
ists; thus it is not easy task to construct the plane-wave 
basis. 

We employ a different form of the vielbein which was 
used at, e.g., Ref . I28l . that is 



B{r)Sl ^^ = 0, 



,3, 



e* — cos 9, e^ = sin 6*, 

el = C{r) sin 6, e^ = C(r) cos ( 



(58) 



The new definition produces the gamma matrices paral- 
lel to the Cartesian coordinate (0:4, 0:5). Non- vanishing 
components of the corresponding spin connections are 



FIG. 5: Spectral flow of the fermion energy of (m, n) — (1,1), 
with some excited energy spectra. The zero crossing spectrum 
has the quantum number K3 = 0, which corresponds to the 
localizing zero mode. 



then 

c^f == 5f,B'ir)cose, Lof = <5^i?'(r) sinfl, 
uf^l-C'{r). (59) 

The Dirac equation is now 

1 9 R' 1 C^' 

i^S^^^d, + z(cos07' + sm9^'){dr + — - -^^) 
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The Dirac gamma matrices should satisfy the anti- 
commutation relations {7,7^} = 2rj and there are 
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FIG. 6: The typical brane solution with (m,n) — (1,3); the profile function and the warp factors (left), and the stress energy 
tensors and the topological charge (right). 



the candidates preserving such Clifford algebra. In most 
of previous studies in 6 dim., their analyses are based 
on the localization on the Abelian vortex and the chi- 
ral representation of gamma matrices. In this represen- 
tation, the spinor can be expanded into the right and 
the left components. The zero mode appears as a eigen- 
state of right or left component. For the massive modes, 
they can be estimated from the mixing of both compo- 
nents. To treat the massive fermionic modes directly, we 
employ the standard representation of the higher dimen- 
sional gamma matrices instead of the chiral one. Since 
the eigenvalues depend on the properties of the back- 
ground brane solutions, zero modes should appear as a 
special case of the massive modes. The standard repre- 
sentation of the gamma matrices in 6 dim. can be defined 
as 
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(61) 



where /„ means the identity matrix of the dimension n. 
The six dimensional spinor '5 can be decomposed into the 
four dimensional and the extra space-time components 

vI/(x^r,0) = v(^^)(^5^[;;;^j) . (62) 

Here the four dimensional part ip{x'^) is the solution of 
the corresponding Dirac equation on the brane 



i^^df^ij} = wTp 
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in which the eigenvalues w indicate an 4 dim. effective 
mass of the fermions. Substituting the ansatz (j62p and 
([55)1 into the Dirac equation (|55p yields the equations for 

- Be-'\dy - ^-J-+Ay)V + Bmf- $U = wU 

Be'\dy+''-J-+Ay)U -BmT-$V ^wV (64) 



where Ay := 2B'/B - (1 - C')/2C. Here we have intro- 
duced the dimensionless coupling constant m := M/eF 
(and uj := uj/eF). In order to eliminate such induced po- 
tential Ay from the equations, it is convenient to replace 
the cigenfunctions {U, V) into (m, v) [23] by 
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The eigenprobleni thus reads 
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H is regarded as a effective Hamiltonian of the 
If background is flat {i.e.,B = 1,5' = 0, C = 
corresponding Hamiltonian easily reads 



H, 



flat 



-7'7'54 - tS'^s 



- 7 rriT ■ i 



where the partial derivatives 84, d^ in Cartesian 
nate are defined as 84 ± id^ — e^^^{dy ± idg/y). 
taining the concise form ([67[) . we introduce and 
additional component of the gamma matrix 7^ := 
One easily confirms that ifflat commutes with 
spin" 

7^ T^ 

K3:=h + — + n—, I3 := 0:4^5 - X5P4 
and also "time-reversal operator" 
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(69) 



where C is the charge conjugation operator. We empha- 
size that these operators also commute with the Hamil- 
tonian in the curved space-time. As a consequence, 
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FIG. 7: Fermion number density for the background skyrmion 
with {m,n) — (1,3) where the coupUng constant m = 0.76. 
The states Kg = 0, ±1 are locaUzed on the brane core. Doubly 
degenerate K3 = ±1 states are more localized (they exhibit 
the non-zero value at the origin). 



FIG. 8: Spectral flow of the fermion energy of (m, n) = (1,3), 
with some excited energy spectra. The zero crossing spectra 
have the quantum numbers K3 — 0, ±1, which correspond to 
the localizing zero modes. 



eigenstates are specified by the magnitude of the grand- 
spin, i.e.. 



K:, 



0,±l,±2,±3---for odd n 
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iiTg = ±-,±-,±-,-- -for even n. (70) 

Since the Hamiltonian is invariant under time reverse, 
the states of zLK^ are degenerate in energy. 

In order to treat the eigenproblem (j65|) . we employ 
the method which was originally proposed by Kahana- 
Ripka [33] for solving the Dirac equation with non-linear 
chiral background. Wc construct the plane-wave basis in 
large circular box with radius _D as a set of eigenstates 
of the flat, unperturbed (B = 1, S' = 0, C = r, / = tt) 
Hamiltonian i.e., Hq = —j^j'^d4 — ^^j^d5 — j^mT^. The 
solutions are 



0,up 



Mk, 



-iVL 



<,. Jp 



i{hy)e 



i(p-i)e 



Ojdown 



<J' 



9+1 



ihy) 



,i{q+l)e 



(71) 



with 



Mk,= 




-1/2 


Nu^ 




-1/2 


and a;+<o,a;->c 
to) . The momen 


= l.<<0'<>0 - - 


-sgn(efc)fc/(efc + 
{) are discretized 



by the boundary conditions 



(72) 



1 71 1 77- 

where p := X3 -|- - - -, g := iiTg - - + -. The orthogo- 
nality of the basis is then satisfied by 



drrJv{kir)Jy(kjr) 



/ drrJ^±i{kir)J^±i{kjr) 
Jo 



D^ 



[J,±iihD)]^ iy:^K3± 



(73) 



Expanding the eigenstates of Eq. (p5|) in terms of the 
plane- wave basis, the eigenproblem reduces to the sym- 
metric matrix diagonalization problem. A special care is 
taken on the estimation of the matrix element from the 
Hamiltonian ([55)1 . In order to hold the Hermiticity of the 
matrix, the following differential rule is imposed 



|5„ 



dydeC{y)^[^^dy^-idyi^^)4'] (74) 



Once the desired eigenfunctions are obtained, angular 
averaged fermion densities on the brane can be estimated 
as follows 

R{y) = N{y)p{y) (75) 

p{y) := / d0[u\y, 0)u{y, 9) + v\y, e)v{y, 9)] (76) 



where 



N{y) = exp 



-41nB(y) 



dy' 

C{y') 



(77) 
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FIG. 9; The brane solutions with large parameter space are plotted. 



For a numerical convenience, we divide the extra space- 
time into two domains. We introduce an effective size 
of brane D* that the stress energy tensors are finite in- 
side. Outside this domain, no brane exists and only the 
warping of the geometry is effective. We assume that the 
fermions we should observe in our 4D space-time are per- 
fectly trapped inside and never leak outside. We set the 
size as D* = D and tentatively choose the value D = 10, 
which is larger than the distribution of the stress energy 
tensors and the topological charge (see Fig[2]). Of course 
outside of the domain is not Minkowski space so that the 
fermions at the boundary feel the effects. The numerical 
results necessarily depend on the choice of D. This crude 
approximation works well if the effects of the geometry 
do not affect to the localized modes of the fermions. 

For the numerical analysis except for the results of 
Figs. |9llD] . we use the brane solutions with following pa- 
rameter set for the background fields 

{■m,n) = (1,1) : 

a = 0.5, Ab = -0.1, il = 0.220525915, 
/'(O) = 0.7292182131078184; (78) 

{m,n) = (1,3) : 

a = 0.11, Afc = -0.1, fL = 0.49588 , 
f^^\0) = 0.1913003760462098. (79) 



The numerical results are shown in Figs. 

Figd] shows the density p{y) for the background 
skyrmion of {m,n) — (1, 1). We present a tower of the 
massive modes as well as the ground state. As one eas- 
ily see that only the lowest mode is peaked on the brane 
while all other modes escape from the core thus we can- 
not observe them. 

In Figdl we plot the densities for some ranges of cou- 
pling constant m. As you see in the figure, only zero 
mode (or quasi zero modes) are localized and the others 
go far away from the brane. In Fig [SI we plot some crucial 
eigenvalues for varying m. The result clearly exhibits so 
called the level crossing picture [37l|. Spectral flow is de- 
fined as the number of eigenvalues of Dirac Hamiltonian 
that cross zero from below minus the number of eigenval- 
ues that cross zero from above for varying the properties 



of the background fields. The spectral flow exactly co- 
incides with the topological charge and the zero modes 
are emerged when they cross the zero. In this case, for 
varying coupling constant m, one zero mode is observed. 

Libanov,Troitsky have discussed relation between the 
topological charge and the fermionic generation [25J . 
Four-dimensional fermions appear as zero-modes trapped 
in the core of the global vortex with winding number 
three. We shall investigate this speculation with our so- 
lution for (m,n) = (1,3). The background brane solu- 
tion is given in Fig[51 As is expected, we observe three 
localized solutions (FigjT]) which can be regarded as the 
generations of the fermions. A main difference from the 
previous studies is that we use non-linear type of soli- 
ton solutions for constructing the branes. They have a 
richer information of the topology and the spectra ex- 
hibit doubly degenerate ground states and a higher ex- 
cited state. On the other hand, the linear type of solitons 
like Abelian vortex have the spectra with double degen- 
eracy only. Fermion generation of our universe comprises 
approximate doubly degenerate states and a very high en- 
ergy state which suggests the underlying topology that 
we employs. 

More generally, we explore the fermion level behavior 
for various brane solutions with large parameter space. 
In Figiniwe present the parameter dependence of the so- 
lutions with (m, n) = (1, 1). By using these, we examine 
the fermion energies and show the results in Fig llOi in 
which we present the level crossing for fixed value of rn, Af, 
(thus the each lines become the function of a,fl). The 
emergence of the zero mode is observed and it strongly 
depends on the property of the branes. 

By using the asymptotics of the warp factors pO|) . we 
easily speculate the behaviour N{y) at zero and infin- 
ity [271 



N{y ^ 0) ^ 2/ 



N{y 



oo 



■ exp(4py 



pe2 



-ePy) . 



(80) 



At the far from the core, the densities of the excited 
modes are more enhanced than the localized mode. 
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FIG. 10; Spectral flow of the fermion energy of {m, n) = 
(1,1) for which m — 1.0, At — —0.1 and different values of 
a, fi. Every dots indicates the fermion energy of the different 
background brane solutions. 



IV. CONCLUSION 

In this article, we have proposed new brane solutions 
in six dimensional space-time and have discussed about 
localization of the fermions on them. The brane have con- 
structed by a baby-skyrmion, a generalization of a 0(3) 
non-linear a model, and realizes regular warped compact- 
ification in 6 dim. anti-de Sitter space-time. The metric 
is non-factorizable with the warp factors either exponen- 
tially diverging or decaying in static solutions. But only 
exponentially decaying warp factors can allow gravity lo- 
calization near the brane in the sense that the 6 dim. 
Planck mass takes a finite value. Such solutions can be 
obtained numerically if we impose suitable boundary con- 
ditions at the distance far from the brane and integrate 
in backward. Forward integration method is unstable as 
a unique set of boundary conditions at the origin does 
not produce an unique solution. 

Once the skyrmion branes were found, we studied the 
fermion localization on the skyrmion branes. The Dirac 
equation in 6 dim. curved space-time was constructed 



in terms of introducing the vielbein and six dimensional 
generalization of the gamma matrices. To treat the mas- 
sive modes directly, we employ the standard representa- 
tion of the gamma matrices. In order to solve the eigen- 
problcm, we introduce the plane wave basis in a large 
circular box with radius D. To study the speculation 
that the charge of the background skyrmions is identi- 
fied as the generation of the fermions of our universe, we 
investigated the fermions with charge Q = 1,3. This con- 
jecture may be confirmed in terms of the fermion spectral 
flow. They exactly coincide with the topological charge 
and practically embody the localized modes. We found 
localized fermion modes corresponding to the topologi- 
cal charge. For (m, n) = (1, 1), we observed the solutions 
with some parameter ranges, which means the existence 
of the massive modes as well as the massless one. For 
(tti, n) — (1, 3), three solutions localized on the brane are 
found. They are the lowest modes of doubly degenerate 
and a excited state, which approximately agree with the 
experimental measurements of fermion masses. On the 
other hand, in the Abelian Higgs model one would expect 
the spectra include only double degeneracy. 

In Ref.[lg|, the authors had an elaborate analysis for 
the massive modes in a 5 dim. anti-de Sitter space-time. 
They set a domain that is defined by the effective poten- 
tial, and then solved the Dirac equation with inside and 
out. By imposing the boundary conditions at the sur- 
face, they obtained a number of massive modes. In our 
case, the setup is more crude; for the numerical reason, 
we tentatively assumed that the fermions are confined 
within a hypothetical domain. If the size goes to infinity, 
there is no effect to the zero mode; the others will be de- 
localized. Therefore the massive modes that we observe 
in Figd] might be an artifact of our numerical set up. 
On the other hand, for (m, n) — (1,3) the lower three 
modes still localize on the vicinity at the origin; thus we 
successfully observe the massive localized mode. 
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